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Abstract 

Magic numbers predicted by a 3-dimensional g-deformed harmonic oscillator with u 9 (3) 
D so g (3) symmetry are compared to experimental data for atomic clusters of alkali metals 
(Li, Na, K, Rb, Cs), noble metals (Cu, Ag, Au), divalent metals (Zn, Cd), and trivalent 
metals (Al, In), as well as to theoretical predictions of jellium models, Woods-Saxon and 
wine bottle potentials, and to the classification scheme using the 3n + I pseudo quantum 
number. In alkali metal clusters and noble metal clusters the 3-dimensional g-deformed 
harmonic oscillator correctly predicts all experimentally observed magic numbers up to 1500 
(which is the expected limit of validity for theories based on the filling of electronic shells), 
while in addition it gives satisfactory results for the magic numbers of clusters of divalent 
metals and trivalent metals, thus indicating that u 9 (3), which is a nonlinear extension 
of the u(3) symmetry of the spherical (3-dimensional isotropic) harmonic oscillator, is a 
good candidate for being the symmetry of systems of several metal clusters. The Taylor 
expansions of angular momentum dependent potentials approximately producing the same 
spectrum as the 3-dimensional g-deformed harmonic oscillator are found to be similar to 
the Taylor expansions of the symmetrized Woods-Saxon and "wine-bottle" symmetrized 
Woods-Saxon potentials, which are known to provide successful fits of the Ekardt potentials. 



PACS numbers: 36.40.-c, 36.40.Cg, 03.65.Fd 
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1. Introduction 

Metal clusters have been recently the subject of many investigations (see [p], 0, |3|, |J for 
relevant reviews). One of the first fascinating findings in their study was the appearance 
of magic numbers, analogous to but different from the magic numbers appearing in the 
shell structure of atomic nuclei 0. Different kinds of metallic clusters [alkali metals (Na 
§ 
01 



Li |1(J 0, K |12[], Rb |L3 



cs 10 m m), n ° bie metais ( Cu mm. a § 

Au Hl6|), divalent metals of the IIB group (Zn, Cd) |19| , trivalent metals of the 
III group (Al. In) |2(J] exhibit different sets of magic numbers. The analogy between the 
magic numbers observed in metal clusters and the magic numbers observed in atomic nuclei 
led to the early description of metal clusters in terms of the Nilsson-Clemenger model 



which is a simplified version of the Nilsson model [22], |23fl of atomic nuclei, in which no spin- 
orbit interaction is included. Further theoretical investigations in terms of the jellium model 
23, 25] demonstrated that the mean field potential in the case of simple metal clusters bears 



great similarities to the Woods-Saxon potential of atomic nuclei, with a slight modification 
of the "wine bottle" type gf| |7|. The Woods-S axon potential itself looks like a harmonic 
oscillator truncated at a certain energy value and flattened at the bottom. It should also be 
recalled that an early schematic explanation of the magic numbers of metallic clusters has 
been given in terms of a scheme intermediate between the level scheme of the 3-dimensional 
harmonic oscillator and the square well M. Again in this case the intermediate potential 
resembles a harmonic oscillator flattened at the bottom. 

On the other hand, modified versions of harmonic oscillators |28, 2^] have been recently 
investigated in the novel mathematical framework of quantum algebras [[30], 31], which are 
nonlinear generalizations of the usual Lie algebras. The spectra of g-deformed oscillators 
increase either less rapidly (for q being a phase factor, i.e. q = e tT with r being real) or 
more rapidly (for q being real, i.e. q = e T with r being real) in comparison to the equidis- 
tant spectrum of the usual harmonic oscillator p2| , while the corresponding (equivalent 
within the limits of perturbation theory or WKB-equivalent) potentials p3l [34]] resemble 



the harmonic oscillator potential, truncated at a certain energy (for q being a phase factor) 
or not (for q being real), the deformation inflicting an overall widening or narrowing of the 
potential, depending on the value of the deformation parameter q. 

Very recently, a g-deformed version of the 3-dimensional harmonic oscillator has been 
constructed f35| , taking advantage of the u g (3) D sOq(3) symmetry [5EI, |3~?||. The spectrum of 
this 3-dimensional g-deformed harmonic oscillator has been found to reproduce very well 
the spectrum of the modified harmonic oscillator introduced by Nilsson [p2] , |23|| , without 
the spin-orbit interaction term. Since the Nilsson model without the spin orbit term is 
essentially the Nilsson-Clemenger model used for the description of metallic clusters [pl |, it 
is worth examining if the 3-dimensional g-deformed harmonic oscillator can reproduce the 
magic numbers of simple metallic clusters and, in the case that this is possible, to determine 
potentials giving the same spectrum as this oscillator and compare them to the symmetrized 
Woods-Saxon and "wine-bottle" symmetrized Woods-Saxon potentials, which successfully 

These are the subjects of the present investigation. 



fit |26, 27] the Ekardt potentials 

It is worth mentioning at this point that an effort has been made to describe the magic 
numbers of metal clusters by a quantum number 3n + / [|] , where n is the number of nodes 
in the solution of the radial Schrodinger equation and I is the angular momentum quantum 
number. This approach has been inspired by the fact that degenerate energy levels in the 
hydrogen atom are characterized by the same value of the quantum number n + 1, due to 
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the so (4) symmetry underlying this system, while degenerate energy levels in the spherical 
harmonic oscillator (i.e. the 3-dimensional isotropic harmonic oscillator) are characterized 
by the same value of the parameter 2n + I, due to the su(3) symmetry underlying this 
system. The 3n + 1 quantum number has been used || to approximate the magic numbers 
of alkali metal clusters with some success, and focusing potentials characterized by this 



degeneracy have been determined [3~S|, but no relevant Lie symmetry could be determined 

mm. 

In Section 2 the 3-dimensional g-deformed harmonic oscillator will be briefly described, 
while in Section 3 the magic numbers provided by this oscillator will be compared with the 
experimental data for Na and Li clusters, as well as with the predictions of other theories 
(jellium model, Woods-Saxon and "wine bottle" potentials, classification scheme using the 
3n + I pseudo quantum number). Additional comparisons of magic numbers predicted by 
the 3-dimensional g-deformed harmonic oscillator to experimental data and to the results 
of other theoretical approaches will be made in Section 4 (for other alkali metal clusters and 
noble metal clusters), Section 5 (for divalent group IIB metal clusters), and Section 6 (for 
trivalent group III metal clusters), while in Section 7 potentials giving approximately the 
same spectrum as the 3-dimensional g-deformed harmonic oscillator will be determined and 
subsequently compared to the symmetrized Woods-Saxon and "wine-bottle" symmetrized 
Woods-Saxon potentials. Finally, Section 8 will contain discussion of the present results 
and plans for further work. 

2. The 3-dimensional g-deformed harmonic oscillator 

The space of the 3-dimensional g-deformed harmonic oscillator consists of the com- 
pletely symmetric irreducible representations of the quantum algebra u g (3). In this space 
a deformed angular momentum algebra, so g (3), can be defined The Hamiltonian of 



the 3-dimensional g-deformed harmonic oscillator is defined so that it satisfies the following 
requirements: 

a) It is an so g (3) scalar, i.e. the energy is simultaneously measurable with the g-deformed 
angular momentum related to the algebra so g (3) and its z-projection. 

b) It conserves the number of bosons, in terms of which the quantum algebras u g (3) and 
so 9 (3) are realized. 

c) In the limit q — > 1 it is in agreement with the Hamiltonian of the usual 3-dimensional 
harmonic oscillator. 

It has been proved that the Hamiltonian of the 3-dimensional g-deformed harmonic 
oscillator satisfying the above requirements takes the form 

H g = tku {[N)q N+1 

where A" is the number operator and is the second order Casimir operator of the algebra 
so g (3), while 

« 2 > 

is the definition of g-numbers and g-operators. 

The energy eigenvalues of the 3-dimensional g-deformed harmonic oscillator are then 




E q (n,l) = |Mg" +1 - q{<1 [2 f 1] m + 1]} . (3) 
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where n is the number of vibrational quanta and / is the eigenvalue of the angular momen- 
tum, obtaining the values I — n, n — 2, . . . , or 1. 

In the limit of q — > 1 one obtains lim q ^iE q (n, I) = hu>on, which coincides with the 
classical result. 

For small values of the deformation parameter r (where q = e T ) one can expand eq. (3) 
in powers of r obtaining |35| 



E q (n, I) = hajQ-n — Tiloqt (1(1 + 1) — n(n + 1)) 

- huj T 2 (l(l + 1) - ^n(n + l)(2n + 1)) + 0(r 3 ). (4) 
The last expression to leading order bears great similarity to the modified harmonic 



oscillator suggested by Nilsson |22|, |23| (with the spin-orbit term omitted) 



V = ^hup 2 - hun'(L 2 - < L 2 > N ), p = r\l^^, (5) 




where 



< V >^ (6) 



The energy eigenvalues of Nilsson's modified harmonic oscillator are [22, 23 1 

E ni = hun - hup! + 1) - ^n(n + 3)^ . (7) 



It has been proved that the spectrum of the 3-dimensional g-deformed harmonic os- 



cillator closely reproduces the spectrum of the modified harmonic oscillator of Nilsson. In 
both cases the effect of the 1(1 + 1) term is to flatten the bottom of the harmonic oscillator 
potential, thus making it to resemble the Woods-Saxon potential. 

The level scheme of the 3-dimensional g-deformed harmonic oscillator (for Huq = 1 and 
t = 0.038) is given in Table 1, up to a certain energy. Each level is characterized by the 
quantum numbers n (number of vibrational quanta) and / (angular momentum). Next 
to each level its energy, the number of particles it can accommodate (which is equal to 
2(2/ + 1)) and the total number of particles up to and including this level are given. If the 
energy difference between two successive levels, which we shall denote by 5, is larger than 
0.39, it is considered as a gap separating two successive shells and the energy difference is 
reported between the two levels. In this way magic numbers can be easily read in the table: 
they are the numbers appearing above the gaps, written in boldface characters. 

Additional level schemes of the 3-dimensional g-deformed harmonic oscillator are given 
in Table 2 (for r = 0.020 and energy gap S = 0.20) and in Table 3 (for r = 0.050 and 
energy gap 5 = 0.38). The following remarks are now in place: 

i) Small magic numbers do not change much as the parameter r is varied (taking positive 
values), while large magic numbers get more influenced by the parameter modification. In 
general, the ordering of the levels does not change rapidly with the value of the parameter 
r (for r > 0). 

ii) Rapid change of the magic numbers as a function of r occurs when r takes negative 
values, but this case is irrelevant to the contents of the present work. 
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iii) Magic numbers are influenced nore drastically by the value of the energy gap S. If in 
the spectrum obtained for a given value of the parameter r the energy separation between 
two successive levels is only slightly smaller than the energy gap 5, this can be considered 
as an indication of the presence of a "secondary" magic number. (See the end of Section 3 
for specific examples.) 

3. Sodium and lithium clusters 

The magic numbers provided by the 3-dimensional g-deformed harmonic oscillator in 
Table 1 are compared to available experimental data for Na clusters 0, 0, § 3 and Li 
clusters fT0| , |TT| in Table 4 (columns 2-7). Some preliminary results concerning Na clusters 
have already been given earlier in Ref. H40| . The following comments apply: 

i) Only magic numbers up to 1500 are reported, since it is known that filling of electronic 
shells is expected to occur only up to this limit ||. For large clusters beyond this point 
it is known that magic numbers can be explained by the completion of icosahedral or 
cuboctahedral shells of atoms |J . 

ii) Up to 600 particles there is consistency among the various experiments and between 
the experimental results in one hand and our findings in the other. 

iii) Beyond 600 particles the results of the four experiments, which report magic numbers 
in this region, are quite different. However, the results of all four experiments are well 
accommodated by the present model. In addition, each magic number predicted by the 
model is supported by at least one experiment. 

In Table 4 the predictions of three simple theoretical models || (non-deformed 3- 
dimensional harmonic oscillator (column 10), square well potential (column 9), rounded 
square well potential (intermediate between the previous two, column 8) ) are also reported 
for comparison. It is clear that the predictions of the non-deformed 3-dimensional harmonic 
oscillator are in agreement with the experimental data only up to magic number 40, while 
the other two models give correctly a few more magic numbers (58, 92, 138), although they 
already fail by predicting magic numbers at 68, 70, 106, 112, 156, which are not observed. 

It should be noticed at this point that the first few magic numbers of alkali clusters 
(up to 92) can be correctly reproduced by the assumption of the formation of shells of 



atoms instead of shells of delocalized electrons |4T| , this assumption being applicable under 
conditions not favoring derealization of the valence electrons of alkali atoms. 

Comparisons among the present results, experimental data for Na clusters (by Martin 
et al. |J (column 2) and Pedersen et al. || (column 3)), experimental data for Li clusters 
(Brechignac et al. [JT0|] (column 4) ), and theoretical predictions more sophisticated than 
these reported in Table 4, can be made in Table 5, where magic numbers predicted by 
various jellium model calculations (columns 5-8, || 0, ||, [42]), Woods-Saxon and wine 
bottle potentials (column 9, |4l|), as well as by a classification scheme using the 3n + I 
pseudo quantum number (column 10, ||) are reported. The following observations can be 
made: 

i) All magic numbers predicted by the 3-dimensional g-deformed harmonic oscillator are 
supported by at least one experiment, with no exception. 

ii) Some of the jellium models, as well as the 3n + I classification scheme, predict magic 
numbers at 186, 540/542, which are not supported by experiment. Some jellium models also 
predict a magic number at 748 or 758, again without support from experiment. The Woods- 
Saxon and wine bottle potentials of Ref. [|3| predict a magic number at 68, for which no 



experimental support exists. The present scheme avoids problems at these numbers. It 
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should be noticed, however, that in the cases of 186 and 542 the energy gap following them 
in the present scheme is 0.329 and 0.325 respectively (see Table 1), i.e. quite close to the 
threshold of 0.39 which we have considered as the minimum energy gap separating different 
shells. One could therefore qualitatively remark that 186 and 542 are "built in" the present 
scheme as "secondary" (not very pronounced) magic numbers. 
4. Other alkali metals and noble metals 

Experimental data for various alkali metal clusters [Li (llOll, column 2), Na (RBfl, column 
3), K (pj, column 4), Rb (pf, column 5), Cs ([0, 0, 
[Cu ([PJ, column 7), Ag ([O] in column 8 and 



column 6)] and noble metal clusters 
in column 9), Au 



nj, column 10)] are 

reported in Table 6, along with the theoretical predictions of the 3-dimensional g-deformed 
harmonic oscillator given in Table 1. The following comments apply: 

|j~3|l , Cu |16|], Ag Au [|16|], what is seen experimentally is cations 
Au 



i) In the cases of Rb 
of the type Rb^, Cu^-, 



Ag^y, Aujy, which contain iV atoms each, but N — 1 electrons. The 
magic numbers reported in Table 6 are electron magic numbers in all cases. 

ii) All alkali metals and noble metals give the same magic numbers, at least within 
the ranges reported in the table. For most of these metals the range of experimentally 
determined magic numbers is rather limited, with Na |J, Cs |7], i~4"| , Li ||10|| , and Ag ]T8| 
being notable exceptions. 

Cs 0, [TJ, Li |TO 



hi) The magic numbers occuring in Na 



and Ag ITH] are almost 



identical, and are described very well by the 3-dimensional g-deformed harmonic oscillator 
of Table 1. The limited data on K, Rb, Cu, Au, also agree with the magic numbers of Table 
1. 

5. Divalent metals of the IIB group 

For these metals the quantities determined experimentally Jl9| are numbers of atoms 
exhibiting "magic" behaviour. Each atom has two valence electrons, therefore the magic 
numbers of electrons are twice the magic numbers of atoms. The magic numbers of electrons 
for Zn and Cd clusters [^] are reported in Table 7 (in columns 4 and 5 respectively), along 
with the magic numbers predicted by the 3-dimensional g-deformed harmonic oscillator for 
two different parameter values (given in Tables 1 and 2 and reported in columns 1 and 2 
respectively), and the magic numbers given by a potential intermediate between the simple 
harmonic oscillator and the square well potential (|19]], column 3). The following comments 
can be made: 

i) The experimental magic numbers for Zn and Cd [[H]] are almost identical. Magic 
numbers reported in parentheses are "secondary" magic numbers, while the magic numbers 
without parentheses are the "main" ones, as indicated in 0. 

ii) In column 1 of Table 7 magic numbers of the 3-dimensional g-deformed harmonic 
oscillator with r = 0.038 and energy gaps larger than 0.26 are reported. Decreasing the 
energy gap 5 considered as separating different shells from 0.39 (used in Table 1) to 0.26 
(used in Table 7) has as a result that the numbers 70 and 106 become magic, in close 
agreement with the experimental data. Similar but even better results are gotten from the 
3-dimensional g-deformed harmonic oscillator of Table 2, reported in column 2 of Table 
7. This oscillator is characterized by r = 0.020, while the energy gap S between different 
shells is set equal to 0.20 . We observe that the second oscillator predicts an additional 
magic number at 112, in agreement with experiment, but otherwise gives the same results 
as the first one. We remark therefore that the general agreement between the results given 
by the 3-dimensional g-deformed harmonic oscillator and the experimental data is not 
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sensitively dependent on the parameter value, but, in contrast, quite different parameter 
values (r = 0.038, r = 0.020) provide quite similar sets of magic numbers (at least in the 
region of relatively small magic numbers). 

hi) Both oscillators reproduce all the "main" magic numbers of Zn and Cd, while the 
intermediate potential between the simple harmonic oscillator and the square well potential, 
reported in column 3, reproduces all the "main" magic numbers except 108. 

6. Trivalent metals of the III group 

Magic numbers of electrons for the trivalent metals Al and In J2U] are reported in Table 



7 (in columns 7 and 8 respectively), along with the predictions of the 3-dimensional q- 
deformed harmonic oscillator of Table 3 (column 6). The following comments can be made: 
i) It is known |], ^] that small magic numbers in clusters of Al and In cannot be 
explained by models based on the filling of electronic shells, because of symmetry breaking 



caused by the ionic lattice ||20|| , while for large magic numbers this problem does not exist. 

ii) The 3-dimensional g-deformed harmonic oscillator of Table 3 provides the magic 
numbers reported in column 6 of Table 7. These magic numbers agree quite well with the 
experimental findings, with an exception in the region of small magic numbers, where the 
model fails to reproduce the magic numbers 164 and 198, predicting only a magic number 
at 186. In addition the oscillator predicts magic numbers at 398, 890, 1074, which are not 
seen in the experiment reported in column 7. 

7. Potentials corresponding to the 3-dimensional g-deformed harmonic os- 
cillator 

As we have seen in the previous sections, the 3-dimensional g-deformed harmonic os- 
cillator describes successfully the magic numbers of several metallic clusters. On the other 
hand, it is known that metallic clusters are successfully described by the Ekardt potentials 
|2"4"|| (for which analytical expressions are lacking), which have been recently parametrized in 
terms of symmetrized Woods-Saxon and "wine-bottle" symmetrized Woods-Saxon poten- 
tials [^, ^7j] (for which analytical expressions are known). Therefore the following questions 
are created: 

a) Is it possible to determine some potentials which, when put into the Schrodinger 
equation, will provide approximately the same spectrum as the 3-dimensional g-deformed 
harmonic oscillator? 

b) If such potentials can be found, how do they compare with the symmetrized Woods- 
Saxon and "wine-bottle" symmetrized Woods-Saxon potentials? 

Question a) is a standard problem of inverse scattering |[44}| . Classical potentials giving 
approximately the same spectrum as the one-dimensional g-deformed harmonic oscillator 



have been determined either through use of standard perturbation theory |3"3 |, or within 
the limits of the WKB approximation fl34 |. 



In what follows we are going to determine potentials giving approximately the same 
spectrum as the 3-dimensional g-deformed harmonic oscillator by using the method of Ref. 



33[| , i.e. perturbation theory. According to this method, a potential of the form 

V = V + kx 2 + Ax 4 + fix 6 + £x 8 + ■ ■ ■ (8) 

corresponds, in first-order perturbation theory and keeping terms up to x 8 only, to a spec- 
trum 

E = e + k + 3A + 15/i + 105£ + (2k + 6A + 40/x + 280£)^ 
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+ (6A + 30/x + + (20// + 140£)™ 3 + 70£n 4 . (9) 

The second term in Eq. (8) corresponds to the usual harmonic oscillator. For appropriate 
values of the numerical coefficients k, A, fi, £, the rest of the terms can be considered as 
perturbations to the harmonic oscillator. 

It is clear that this method can be applied in cases in which the spectrum under study 
depends on only one quantum number, the number of excitation quanta n. In the case of 
the 3-dimensional g-deformed harmonic oscillator (Eq. (3)), though, the spectrum depends 
on an additional quantum number, the angular momentum I. One way out is to determine 



a /-dependent equivalent potential, as it is done in several branches of physics [|45|, |46[| . In 
order to do this, for each possible value of I (I = n, n — 2, n — 4, . . . 1 or (see Eq. (3)) 
one determines the energy as a function of n only and then calculates the corresponding 
potential. 

In the case of the 3-dimensional g-deformed harmonic oscillator the energy spectrum 
E q (n,l) for the various possible values of the angular momentum I (I = n, n — 2, n — 4, 
. . . 1 or 0) can be put in the form: 





E q (n,n) = hu [n\ q 2 for I = n, 




(10) 




E q (n, n-2) = fnu (q 2 [n - l] g2 + q 2lA j for I = n - 2, 




(11) 


E q ( 


n, n - 4) = huj (q 4 [n - 2} q 2 + q 2{n '^ [2]^) for I = n - 


4, 


(12) 


E q ( 


n, n- 6) = hco {q 6 [n - 3} q 2 + q 2 ^- 2 \3} q 2^j for l = n- 


6, 


(13) 




E q (n, 3) = hu (q 8 [[n - 3)) q2 + q~ A [[5]] q 2 - q- 2 [[3]] q 2 + l) 








= E q (n,0) -hu (q 6 - l)(l + g~ 4 ) for / = 3, 




(14) 


E q (n, 2) = hcu (q%n - 2)) q2 + g~ 2 [[3]] 9 2 - l) = E q (n, 0) - huj (q 4 - q~ 2 ) 


for 


I = 2, 




) = huo (q 4 [[n - !}} q 2 + l) = E q {n, 0) - hco (q 2 - 1) for 




(15) 


E q (n,l 


1 = 1, 


(16) 




E q (n, 0) = ^ti; g 2 [[n]] g 2 for I = 0, 




(17) 



where by [n] q are denoted the g-numbers of Eq. (2), which are symmetric under the exchange 
q <-> g -1 , while by [[n]] q are denoted the g-numbers 

[ML = (is) 

which are not symmetric under the exchange q <-> g _1 . For all of these equations it is clear 
that they reduce to the classical expression E(n) = hu^n in the limit g — » 1. 

We then consider the Taylor expansions for these energy expressions. By comparing 
them to Eq. (9) and equating the coefficients of the various powers of n (up to n 4 ) we 
determine in each case the coefficients k, X, \i, £. Substituting these coefficients in Eq. 
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(8) we determine for each case the corresponding potential, keeping terms up to r 4 (where 
q = e T ). The first few cases are: 

^n, _ ( l 4r 4 \ (l r 2 4r 4 \ , 

v 2 2 15 / V 2 2 15 y 

-(t-it) i4+ (s-S) 16 for <="• (19 > 



V(x)i =n _ 2 1 A 7r 2 8r 3 8r 4 \ /l n 3r 2 10r 3 44r 4 , 

V =- - + 4r + — + — + — + - + 2r+ — + — - + — — x 
huj \2 2 3 5 / \2 2 3 15 / 



2 



1 , r r 3 r 4 \ / 1 r 2 4r 4 \ /r 2r 3 r 4 
V 2 2 6 15 i + l2 2 + 15 i + 1 6 9 6 J 
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+ ( — - — ) x 6 + ( — + — | x 8 for / = 0, (21) 
V 30 45 / V 210 210 / 

v^(x), =1 = \/(x) i=0 - ^ (g 2 - i) 

/ 4r 3 2r 4 \ 

~V(x) l=0 -hu l2r + 2r 2 + — + — J for 1 = 1, (22) 

V(x) l=2 = V(x) l=0 -huj (q 4 -q' 2 ) 
~ \/(a;) /=0 - 7kj (6t + 6r 2 + 12r 3 + 10r 4 ) for Z = 2, (23) 

V(x) z=3 = V(x) l=0 - ku (q e - 1)(1 + g" 4 ) 
^ V(x) i=0 - ^ (12r + 12r 2 + 48r 3 + 44r 4 ) for / = 3, (24) 



We remark that for small values of r, like the ones occuring in the previous sections, 
the potentials occuring for I = n and I = n — 2 are of the form 

V(x) = V + ax 2 - bx 4 + cx\ (25) 

with a,b,c> 0, while the potentials occuring for I — 0, 1, 2, 3 are of the form 

V{x) = V + ax 2 + bx 4 + cx 6 + dx 8 , (26) 

with a, b,c,d > 0. 

It is instructive at this point to compare these potentials with the symmetrized Woods- 
Saxon potential 

sinh(_R/a) , . 

V S ws(r) = -V K 1 ! < r < oo, 27 

cosh(r/a) + cosh(it/a) 
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and the "wine-bottle" symmetrized Woods-Saxon potential 



/ \ T . ( wr 2 \ smh(R/a) „ 

V WB (r) = -vM + w ) coMr/a) lii Wa y < r < oo, (28) 



which have been used for parametrizing the Ekardt potentials p^| . In order to 

facilitate the comparisons, we consider the Taylor expansions of these potentials 

V S ws(r) 1 1 r 2 



V sinh^/a) 1 + cosh(i?/a) 2(1 + cosh(i?/a)) 2 a 2 
5-cosh(i2/a) r 4 (cosh^/a)) 2 - 28 cosh(i?/a) + 61 r 6 

H 7on/, , ZZUTTTT-vTI Ik ' I 29 ) 



24(1 + cosh( J R/a)) 3 a 4 720(1 + cosh^/a)) 4 

VV s (r) 1 / 1 id a 2 \r 2 



V sinh(i?/a) 1 + cosh( J R/a) \2(1 + cosh(i?/a)) 2 1 + cosh( J R/a) i? 2 / a 2 

/ 5 — cosh(_R/a) w a 2 \ r 4 

~ V24(l + cosh( J R/a)) 3 ~ 2(1 + cosh(i?/a)) 2 i? 2 J ^ 4 

/ (cosh(i?/a)) 2 - 28 cosh(i?/a) + 61 w(5 - cosh(i?/a)) a 2 \ r 6 
+ ^ 720(1 + cosh(i?/a)) 4 24(l + cosh( J R/a)) 3 !R 2 J ^' ' '° ' 

The following comments are now in place: 

i) The Taylor expansions of the symmetrized Woods-Saxon and the "wine-bottle" sym- 
metrized Woods-Saxon potentials, which have been used for fitting the Ekardt potentials 
used for the description of metallic clusters, have the same form as the potentials corre- 
sponding to the 3-dimensional g-deformed harmonic oscillator, i.e. they contain all the even 
powers of the relevant variable (and no odd powers). It is therefore not surprising that the 
3-dimensional g-deformed harmonic oscillator gives a good description of the magic numbers 
of metallic clusters. 

ii) The potentials obtained through the use of perturbation theory are valid near the 
origin (x = 0) and for relatively low values of n. They do not give information about the 
shape of the potential near its edges, or for very large values of n. The determination of 
potentials which will be accurate near their edges remains an open problem. It should 
also be examined if these potentials possess any deeper relation to the quantum algebraic 
symmetry characterizing the 3-dimensional g-deformed harmonic oscillator. For example, 
one could check if these potentials are related to the generators of the relevant quantum 
algebra. The existence of such a relation also remains an open problem. 

iii) For very large values of n, the spectrum gets an exponential form. For example, Eq. 
(10) becomes (for r > 0) 

g2rn g— 2rn g 2Tn 

E q (n,n) = hu e2r _ e _ 2r - ^o e2T _ e _ 2T - (31) 

Potentials with exponential spectra have been considered in Ref. P7| , but also in this case 
only the form of the potential near the origin could be determined. 
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iv) Focusing potentials leading to 3n + I degeneracy of the energy levels (which has 
been found to describe reasonably well the magic numbers of alkali clusters |J) have been 



determined in Ref. [38 . They have the form 



Both of them are of the form 



4 

X 



V(x) = E - AwTJ¥ , (34) 

which corresponds to a Taylor expansion of the form 

V(x) = E-A{x 4 -2x w + 3x w + ■••)• (35) 

We remark that this Taylor expansion bears no similarity to the Taylor expansions of the 
symmetrized Woods-Saxon and "wine-bottle" symmetrized Woods-Saxon potentials, since 
it contains only some of the even powers of the relevant variable and not all of them. 
Indeed, these focusing potentials are known to exhibit the "wine-bottle" feature strongly 
exaggerated |5Bfl , lacking in parallel the flat bottom characterizing the Woods-Saxon and 



Ekardt potentials. The potential U%(r) has, however, the major advantage that it reproduces 
quite well the edge behaviour of the Ekardt potentials ]38 . 
8. Discussion 

The following general remarks can now be made: 

i) From the results reported above it is quite clear that the 3-dimensional g-deformed 
harmonic oscillator describes very well the magic numbers of alkali metal clusters and noble 
metal clusters in all regions, using only one free parameter (q = e T with r = 0.038). It 
also provides an accurate description of the "main" magic numbers of clusters of divalent 
group IIB metals, either with the same parameter value (r = 0.038) or with a different one 
(r = 0.020). In addition it gives a satisfactory description of the magic numbers of clusters 
of trivalent group III metals with a different parameter value (r = 0.050). 

ii) It is quite remarkable that the 3-dimensional g-deformed harmonic oscillator repro- 
duces long sequences of magic numbers (Na, Cs, Li, Ag) at least as accurately as other, 
more sophisticated, models by using only one free parameter (q = e T ). (It should not be 
forgotten at this point that these other models have deep physical roots, while the present 
approach is based on symmetry arguments, which are justified a posteriori by their success- 
ful predictions.) Once the parameter is fixed, the whole spectrum is fixed and no further 
manipulations can be made, the choice of the energy gap 5 being the only exception. How- 
ever, the choice of the value of the energy gap 6 does not influence the order of the energy 
levels, but it is just deciding which energy separations will be considered as corresponding 
to main magic numbers and which not. The successful prediction of the magic numbers 
can be considered as evidence that the 3-dimensional g-deformed harmonic oscillator owns 
a symmetry (the u g (3) D so 9 (3) symmetry) appropriate for the description of the physical 
systems under study. 
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iii) As we have already mentioned, it has been remarked |J that if n is the number of 
nodes in the solution of the radial Schrodinger equation and I is the angular momentum 
quantum number, then the degeneracy of energy levels of the hydrogen atom characterized 
by the same n + I is due to the so(4) symmetry of this system, while the degeneracy 
of energy levels of the spherical harmonic oscillator (i.e. of the 3-dimensional isotropic 
harmonic oscillator) characterized by the same 2n + I is due to the su(3) symmetry of 
this system. 3n + I has been used || to approximate the magic numbers of alkali metal 
clusters with some success, and focusing potentials characterized by this degeneracy have 



been determined [38], but no relevant Lie symmetry could be determined [38, 32]. In view of 
the present findings the lack of Lie symmetry related to 3n + 1 is quite clear: the symmetry 
of the system appears to be a quantum algebraic symmetry (u 9 (3)), which is a nonlinear 
extension of the Lie symmetry u(3). 

iv) The ability of the 3-dimensional g-deformed harmonic oscillator to reproduce cor- 
rectly the magic numbers of several metal clusters is not coming as a surprise, if one con- 
siders potentials giving approximately (within the limits of perturbation theory) the same 
spectrum as this oscillator. The Taylor expansions of these potentials have the same form 
as the Taylor expansions of the symmetrized Woods-Saxon and "wine-bottle" symmetrized 
Woods-Saxon potentials, which successfully fit ^7) the Ekardt potentials |23J, which 
characterize the structure of metal clusters. 

In summary, we have shown that the 3-dimensional g-deformed harmonic oscillator with 
u g (3) D sOq(3) symmetry correctly predicts all experimentally observed magic numbers 
of alkali metal clusters and of noble metal clusters up to 1500, which is the expected 
limit of validity for theories based on the filling of electronic shells. In addition it gives 
a good description of the "main" magic numbers of group IIB (divalent) metal clusters, 
as well as a satisfactory description of group III (trivalent) metal clusters. This indicates 
that u g (3), which is a nonlinear deformation of the u(3) symmetry of the spherical (3- 
dimensional isotropic) harmonic oscillator, is a good candidate for being the symmetry of 
systems of several metal clusters. Furthermore, the Taylor expansions of potentials giving 
approximately the same spectrum as the 3-dimensional g-deformed harmonic oscillator are 
found to have the same form as the Taylor expansions of the symmetrized Woods-Saxon 
and "wine-bottle" symmetrized Woods-Saxon potentials, which sucsessfully fit the Ekardt 
potentials underlying the structure of metal clusters. Naturally, these Taylor expansions 
are valid near the origin. The determination of potentials which will be accurate near their 
edges remains an open problem. It is also an open problem the existence of any deeper 
relation between these potentials and the quantum algebra characterizing the 3-dimensional 
g-deformed harmonic oscillator, as, for example, some relation between these potentials and 
the generators of the quantum algebra. 
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Table 1: Energy spectrum, E q (n,l), of the 3- dimensional g-deformed harmonic oscillator 
(Eq. (3)), for Tiujq = 1 and q = e T with r = 0.038. Each level is characterized by n (the 
number of vibrational quanta) and I (the angular momentum). 2(2/ + 1) represents the 
number of particles each level can accommodate, while under "total" the total number 
of particles up to and including this level is given. Magic numbers, reported in boldface, 
correspond to energy gaps larger than 5 = 0.39, reported between the relevant couples of 
energy levels. 
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Table 1: 



n 


/ 


E q (n,l) 


2(2/ + 1) 


total 


n 


/ 


E q (n,l) 


2(2/+ 1) 
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20 
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6 
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26 
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3.420 


6 
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0.638 






10 
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18 
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4 
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18 
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0.559 
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15.206 


54 
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10 


68 


10 


2 


15.316 


10 
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4 





4.854 


2 


70 
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15.554 


2 
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5 


5 


5.116 


22 


92 


11 


7 
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30 
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12 


10 
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42 
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5 


3 
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14 
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6 
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6.204 


26 
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11 
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22 
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5 


1 
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6 
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14 
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58 
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18 
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11 
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17.385 


14 


926 


7 


7 
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30 
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12 
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17.410 


34 
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6 
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10 
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13 
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46 
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11 
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6 
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14 
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12 
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18 


1206 


9 
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38 
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26 


338 


12 
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10 
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8 
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18 
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14 
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Table 2: Same as Table 1, but with Huq = 1 and q = e T with r = 0.020. The energy gap 
separating different shells has been taken to be 5 = 0.20 . 
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£,(ra,Z) 
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n 


I 


£,(ra,Z) 
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22 
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Table 3: Same as Table 1, but with hu) = 1 and q = e T with r = 0.050. The energy gap 
separating different shells has been taken to be 5 = 0.38 . 
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Table 4: Magic numbers provided by the 3-dimensional (/-deformed harmonic oscillator 
(Table 1), reported in column 1, are compared to the experimental data for Na clusters by 
Martin et al. [6] (column 2), Bj0rnholm et al. [7] (column 3), Knight et al. [8] (column 
4), and Pedersen et al. [9] (column 5), as well as to the experimental data for Li clusters 
by Brechignac et al. ([10] in column 6, [11] in column 7). The magic numbers provided 
[5] by the (non-deformed) 3-dimensional harmonic oscillator (column 10), the square well 
potential (column 9) and a rounded square well potential intermediate between the previous 
two (column 8) are also shown for comparison. See text for discussion. 
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Table 5: Magic numbers provided by the 3-dimensional g-deformed harmonic oscillator 
(Table 1), reported in column 1, are compared to the experimental data for Na clusters 
by Martin et al. [6] (column 2), and Pedersen et al. [9] (column 3), as well as to the ex- 
perimental data for Li clusters by Brechignac et al. [10] (column 4) and to the theoretical 
predictions of various jellium model calculations reported by Martin et al. [6] (column 5), 
Bj0rnholm et al. [7] (column 6), Brack [2] (column 7), Bulgac and Lewenkopf [42] (column 
8), the theoretical predictions of Woods-Saxon and "wine bottle" potentials reported by 
Nishioka et al. [43] (column 9), as well as to the magic numbers predicted by the classi- 
fication scheme using the 3n + I pseudo quantum number, reported by Martin et al. [6] 
(column 10). See text for discussion. 



1 1 

th. 


exp. 


exp. 


exp. 


th. 


th. 


th. 


th. 


th. 


th. 


present 


AT 

Na 


AT 

Na 


T * 

Li 


1 1 

jell. 


1 1 

jell. 


jell. 


1 1 

jell. 


WS 


3n + I 


Tab. 1 


Ref . [6] 


T~» £ [r\~\ 

Ref. [9] 


t~> £ [-i r\~\ 

Ref. [10] 


Ref. [6] 


Ref. [7] 


Ref. [2] 


Ref. [42] 


Ref. [43] 


Ref. [6] 


2 


2 






2 


2 


2 




2 


2 


8 


8 






8 


8 


8 




8 


8 


(18) 


18 






18 


18 








18 


20 


20 






(20) 


20 


20 




20 




34 


34 






34 


34 


34 


34 




34 


40 


40 


40 




(40) 


40 






40 




58 


58 


58 




58 


58 


58 


58 


58 


58 


















68 




92 


90,92 


92 


93 


92 


92 


92 


92 


92 


90 


138 


138 


138 


134 


134 


138 


138 


138 


138 


132 










186 


186 


186 


186 




186 


198 


198±2 


198 


191 


(196) 


196 






198 




254 








254 


254 


254 


254 


254 


252 


268 


263±5 


264 


262 


(268) 








268 




QQfi 
OOO 


o4l±0 


Q/l A 
o44 


Q/1 1 


ooo(oi>o ) 


QQfi 
OOO 


QQfi 
OOO 


QQfi 
OOO 


QQfi 
OOO 


QQO 
OOZ 


440 


443±5 


442 


442 


440 


440 


438,440 


440 


440 


428 














542 


542 




540 


556 


557±5 


554 


552 


562 


556 


556 


556 


562 




676 




680 






676 


676 


676 




670 


694 


700±15 




695 


704 








694 
















758 


748 






832 


840±15 


800 


822 


852 


832 


832 


832 


832 


820 


912 






902 






912 


912 






1012 


1040±20 


970 


1025 






1074 


1074 


1012 


990 


1100 




1120 








1100 


1100 


1100 




1206 


1220±20 














1216 


1182 


1284 






1297 






1284 


1284 






1314 




1310 












1314 




1410 


1430±20 
















1398 


1502 




1500 








1502 


1502 


1516 





21 



Table 6: Magic numbers provided by the 3-dimensional g-deformed harmonic oscillator 
(Table 1), reported in column 1, are compared to the experimental data for clusters of Li 
[10] (column 2), Na [6] (column 3), K [12] (column 4), Rb [13] (column 5), Cs [7,14] (column 
6), Cu [16] (column 7), Ag ([18] in column 8, [16] in column 9), and Au [16] (column 10). 
See text for discussion. 
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Table 7: Magic numbers provided by the 3-dimensional g-deformed harmonic oscillator of 
Table 1 with energy gap 5 = 0.26 (column 1) and of Table 2 (column 2), are compared 
to the experimental data for Zn clusters [19] (column 4) and Cd clusters [19] (column 5), 
as well as to the theoretical predictions of a potential intermediate between the simple 
harmonic oscillator and the square well potential [19] (column 3). In addition, the magic 
numbers provided by the 3-dimensional g-deformed harmonic oscillator of Table 3 (reported 
in column 6) are compared to the experimental data for Al [20] (column 7) and In [20] 
(column 8). See text for discussion. 
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